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At high temperatures, strongly interacting matter becomes a plasma of deconfined 
quarks and gluons. In statistical QCD, deconfinement and the properties of the resulting 
quark-gluon plasma can be investigated by studying the in-medium behaviour of heavy 
quark bound states. In high energy nuclear interactions, quarkonia probe different aspects 
of the medium formed in the collision. So, we first reviewed the fate of quarkonia in the 
different stages of the (dynamical) system produced at the collision. We have then 
presented our present work on the dissociation of the heavy quarkonium states in a hot 
QCD medium by investigating the medium modifications to heavy quark potential. In 
contrast to the usual screening picture, interestingly our theory gives rise the screening 
of the charge, not the range of the potential. 



Introduction: The study of the funda- 
mental forces between quarks and gluons is 
an essential key to the understanding of QCD 
and the occurrence of different phases which 
are expected to show up when going from low 
to high temperatures (T) and/or baryon num- 
ber densities. For instance, at small or vanish- 
ing temperatures quarks and gluons get con- 
fined by the strong force while at high tem- 
peratures asymptotic freedom suggests a quite 
different QCD medium consisting of rather 
weakly coupled deconfined quarks and glu- 
ons, the so-called quark gluon plasma (QGP). 
The anomalous suppression of the J/ip pro- 
duction in heavy ion collisions which has been 
experimentally observed [l[ in the depletion 
of the dilepton multiplicity in the region of 
invariant mass corresponding to the J/ip me- 
son was proposed long time ago as a possi- 
bly unambiguous signal of the onset of decon- 
finement. Matsui and Satz 0] argued that 
charmonium states produced before the for- 
mation of a thcrmalizcd QGP would tend to 
melt in their path through the deconfined 
medium, since the binding (colour) Coulomb 
potential is screened by the large number of 
colour charges. This, in turn, would produce 
an anomalous (with respect to normal nuclear 
absorption) drop in the J/ip yields. 

In this picture it is implicitly assumed that, 
once the charmonium dissociates, the heavy 
quarks hadronizc by combining with light 



quarks only (recombination leading to a sec- 
ondary J/ip production is neglected). This 
assumption is certainly justified at the SPS 
conditions, due to the very small number of 
cc pairs produced per collision (N c5 ~ 0.2 in 
a central collision), but at RHIC (N cS ~ 10) 
and LHC (iV c£ ~ 200) energies it is no longer 
warranted [|[. 

Moreover in a hadronic collisions only 
about 60% of the observed J/ip's are directly 
produced, the remaining stemming from the 
decays of excited charmonium states (notably 
the Xc and the ip'). Since each cc bound 
state dissociates at a different temperature, 
a model of sequential suppression was de- 
veloped, with the aim of reproducing the 
J/ip suppression pattern as a function of 
the energy density reached in the heavy ion 
collision. SPS experimental data for Pb-Pb 
collisions at different ccntralities seem indeed 
to support the dissociation pattern predicted 
by this model. 

The heavy quark pair leading to the J/ip 
mesons are produced in nucleus-nucleus col- 
lisions on a very short time-scale ~ l/2m c , 
where m c is the mass of the charm quark. The 
pair develops into the physical resonance over 
a formation time and traverses the plasma 
and (later) the hadronic matter before leaving 
the interacting system to decay (into a dilep- 
ton) to be detected. This long 'trek' inside 
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the interacting system is fairly 'hazardous' for 
the J /if). Even before the resonance is formed 
it may be absorbed by the nucleons stream- 
ing past it [3|. By the time the resonance is 
formed, the screening of the colour forces in 
the plasma may be sufficient to inhibit a bind- 
ing of the cc Q ■ Or an energetic gluon [f| or 
a comoving hadron could dissociate the reso- 
nance^). 

Quarkonia at finite temperature are an im- 
portant tool for the study of quark-gluon 
plasma formation in heavy ion collisions. 
Many efforts have been devoted to determine 
the dissociation temperatures of QQ states in 
the deconfined medium, using either lattice 
calculations of quarkonium spectral functions 
or non-relativistic calculations based upon 
some effective (screened) potential. 

Lattice studies are directly based on quan- 
tum chromodynamics and should provide, in 
principle, a definite answer to the problem. 
However, in lattice studies the spectral func- 
tions have to be extracted — using rather lim- 
ited sets of data — from the Euclidean (imagi- 
nary time) correlators, which arc directly mea- 
sured on the lattice. This fact, together with 
the intrinsic technical difficulties of lattice cal- 
culations, somehow limits the reliability of the 
results obtained so far, and also their scope, 
which in fact is essentially limited to the mass 
of the ground state in each QQ channel. Po- 
tential models, on the other hand, provide a 
simple and intuitive framework for the study 
of quarkonium properties at finite tempera- 
ture, allowing one to calculate quantities that 
are beyond the present possibilities for lattice 
studies. The main problem of the latter ap- 
proach is the determination of the effective po- 
tential: although at zero temperature the use 
of effective potentials and their connection to 
the underlying field theory is well established, 
at finite T the issue is still open. 

Calculations of the cc and bb dissociation 
temperatures, using different potential mod- 
els based upon the lattice free and internal 
energies, have found on the whole a reason- 
able agreement with the results from the lat- 
tice studies On the other hand, calcula- 
tions of Euclidean correlators using a variety 



of potential models were not able to reproduce 
the temperature dependence of the lattice cor- 
relators. 

A precise quantitative agreement with the 
lattice correlators should not be expected, be- 
cause of uncertainties coming from a variety 
of sources. Not only the determination of the 
effective potential is still an open question but 
there are also issues tied, e.g., to relativistic ef- 
fects, to the thermal width of the states or to 
the contribution of radiative corrections. On 
the other hand, lattice correlators are also af- 
fected by their own uncertainties. These may 
be due to the use of different lattices (isotropic 
or anisotropic); to the finite size of the box, 
which might significantly alter the continuum 
part of the spectrum, although calculations 
with boxes of different sizes show discrepan- 
cies below 1% [5,6]; or to artifacts in the con- 
tinuum region of the spectral functions due to 
the finite lattice spacing [5]. 

Recently, Umeda and Alberico (nj EH have 
shown that the lattice calculations of meson 
correlators at finite temperature contain a 
constant contribution, due to the presence of 
zero modes in the spectral functions. These 
contributions cure most of the previously ob- 
served discrepancies with lattice calculations, 
supporting the use of potential models at finite 
temperature as an important tool to comple- 
ment lattice studies. 

Actually, even if the potential supports the 
existence of bound states, other physical pro- 
cesses may lead to the dissociation of the 
quarkonium. First, if the QQ binding energy 
is lower than the temperature — and assum- 
ing that the quarkonia have reached the ther- 
mal equilibrium with the plasma — a certain 
fraction of their total number will be ther- 
mally excited to resonant states according to 
a Bosc-Einstcin distribution: such a process is 
referred to as thermal dissociation. Further- 
more, the collisions with the gluons and the 
light quarks of the plasma may lead to the 
collisional dissociation of the quarkonium Q . 

Binoy and Menon revisited the J/ip sup- 
pression due to gluonic bombardment in an 
expan ding quark-gluon plasma in a series of 
works [12j ■ First they neatly incorporated the 
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crucial effects arising from gluon fugacity, rel- 
ative g~ ip flux, and J/^meson formation time 
and then used these effects in the formulation 
of the gluon number density, velocity- weighted 
cross section, and the survival probability in 
an equilibrated static QGP. This formulation 
have been used to study the pattern of J/ijj 
suppression in the central rapidity region at 
RHIC/LHC energies. Later they explicitly 
take into account the effect of hydrodynamic 
(longitudinal) expansion profile on the gluonic 
breakup of J/^'s in an (chemically) equilibrat- 
ing expanding QGP. A novel type of partial- 
wave interference mechanism is found to op- 
erate in the modified dissociation rate. Fi- 
nally, this formulation has been applied to the 
case when the medium is undergoing cylindri- 
cally symmetric transverse expansion. Com- 
pared to the case of longitudinal expansion 
the new graph of survival probability develops 
a rich structure at RHIC, due to a competi- 
tion between the transverse catch-up time and 
plasma lifetime. 

Of course, if the process g + J/ip — ► c + c, 
discussed above, can lead to the dissociation 
of the charmonium, the same reaction can also 
occur in the opposite direction. Hence a con- 
sistent calculation of J/-0multiplicity implies 
the solution of a kinetic rate equation inte- 
grated over the lifetime of the QGP phase in 
which both processes (dissociation and recom- 
bination) enter 0] This is of relevance be- 
cause, as mentioned above, the usual assump- 
tion in considering the J/?/>supprcssion as a 
signature of deconfinement is that its produc- 
tion can occur only in the very initial stage 
of the collision. Really, if at SPS the role 
played by recombination is numerically negli- 
gible, this is no longer true at RHIC as pointed 
out in Ref. ||. 

However, there is a generic question, quite 
often asked, is that whether we can distinguish 
between the two mechanisms of dissociation 
(colour screening and collision with hard glu- 
ons) mainly operating in the deconfined phase. 
Binoy and Srivastava [l3[ have shown that 
while the gluonic dissociation of the J/ ipis al- 
ways possible, the Debye screening is not ef- 
fective in the case of small systems at RHIC 



energies. For the larger systems, the Debye 
screening is more effective for lower transverse 
momenta, while the gluonic dissociation dom- 
inates for larger transverse momenta. At LHC 
energies the Debye screening is the dominant 
mechanism of J / ip suppression for all the cases 
and momenta studied. As an interesting re- 
sult, they found the gluonic dissociation to be 
substantial but the Debye screening to be in- 
effective for T suppression at the LHC energy. 

In the context of J/f/'supprcssion, Langcvin 
dynamics seems to be almost as important 
as other mechanisms invoked to explain the 
RHIC and LHC data. Binoy and Menon Q 
considered the Brownian motion of a cc pair 
produced in the very early stage of a QGP. 
They found that, in the weak coupling regime, 
both the time scales associated with the posi- 
tional swelling (t x ) and approach to ionization 
(te) arc positive and less than the frictional 
relaxation time (7~ 1 ). Hence Brownian move- 
ment can cause a genuine break-up of the cc 
bound by swelling it substantially or Langevin 
dynamics can cause cc to ionize after a time 
span te- On the other hand, in the strong 
coupling case, t x is imaginary and te are neg- 
ative, i.e., unphysical. Hence random force 
plus diffusion cannot cause the cc bound state 
to dissociate. 

While the short and intermediate distance 
(rT < 1) properties of the heavy quark in- 
teraction is important for the understanding 
of in-medium modifications of heavy quark 
bound states, the large distance property 
of the heavy quark interaction which is 
important for our understanding of the bulk 
properties of the QCD plasma phase viz. 
the equation of state. In all of these studies 
deviations from perturbative calculations and 
the ideal gas behaviour are expected and were 
indeed found at temperatures which are only 
moderately larger than the deconfinement 
temperature. This calls for quantitative 
non-pcrturbative calculations. The phase 
transition in full QCD will appear as an 
crossover rather than a 'true' phase transition 
with related singularities in thermodynamic 
observables (in the high-temperature and 
low density regime) a cross-over, it can be 
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reasonable to assume that the string-tension 
does not vanish abruptly above T c . So we [l5[ 
decide to investigate in our present work 
what happens to the different quarkonium 
states if one corrects with a dielectric func- 
tion encoding the effects of the deconfined 
medium, the full Cornell potential and not 
only its Coulomb part as usually done in the 
literature. We found that with this choice 
medium effects give rise to a long-range 
Coulomb potential with a reduced effective 
charge (inversely proportional to the square 
of the Debye mass) of the heavy quark, at 
variance with its usual Debye-screened form 
employed in most of the literature. With such 
an effective potential we investigate the ef- 
fects of different possible choices of the Dcbyc 
mass on the dissociation temperature of the 
different quarkonium states. Since a Coulomb 
interaction always admits bound states, a 
criterion has to be adopted to define such 
a dissociation temperature: a state is then 
considered "melted" when its binding energy 
becomes of the same order as the temperature. 

The Debye masses in hot QCD: 

Adequate knowledge of Debye mass is in- 
deed needed to study the medium modifica- 
tions to heavy quark potential. The Debye 
mass in QCD unlike QED is generically non- 
perturbative and gauge invariant [16]. The De- 
bye mass at high temperature in the leading- 
order in QCD coupling is known from long 
time and is perturbative [ItJ ■ The Debye mass 
in leading-order from the polarization tensor 
of a gauge boson derived from the HTL ap- 
proach can also be obtained from the trans- 
port theory (l8j . One cannot naively gen- 
eralize the definition of Debye mass in QED 
to QCD due to the non-abelian nature of the 
theory. Rebhan [l9[ has defined Debye mass 
through the relevant pole of the static quark 
propagator instead of the zero momentum 
limit of the time-time component of the gluon 
self-energy. The Debye mass thus defined 
comes out to be gauge independent follows 
from the fact that the pole of the self-energy is 
independent of choice of gauge. Braaten and 
Nieto [2(j computed the Debye screening mass 



for QGP at high temperature to the next-to- 
lcading-ordcr in QCD coupling from the cor- 
relator of two Polyakov loops which agreed to 
the HTL result of [lj]. 

Arnold and Yaffe [lj| pointed out that the 
contribution of order (g 2 T) to the Debye mass 
in QCD needs the knowledge of the non- 
perturbative physics of confinement of mag- 
netic charges and a perturbative definition of 
the Debye mass as a pole of gluon propaga- 
tor no longer holds. They showed how one 
can define Debye mass in QCD in a mani- 
festly gauge invariant manner (in vector-like 
gauge theories with zero chemical potential). 
In the work of Kajantic et. al [2l| . the 
non-perturbative contributions of 0(g 2 T) and 
0(g 3 T) have been determined from 3-D ef- 
fective field theory which we consider in the 
present work. At high temperatures and zero 
chemical potential Debye mass can be ex- 
panded in a power series in QCD coupling 21]: 



m D 



LO 



Ng 2 T 

47T 



111 



g 2 T 



,g 2 T + d N Nf g 3 T + 0(g 4 T) , (1) 



where is the leading-order result [l7j ]. 

The coefficient d,N,N f have the following de- 
pendence on the number of colors N and fla- 
vors Nf as: 



>N 



IN.N, 



y/N/3 + N f /6 



(2) 



where the values of c N , &w have been obtained 
by fitting the results with the physical 4D fi- 
nite temperature QCD [2lj : 



SU(3) 



2.46 ±0.15 b 



N 



-0.49 ±0.15(3) 



The number c N captures the non-perturbative 
3-D effects, while the d,N,N f is related to the 
choice of scale in m^°. We employ the two- 
loop expression for the QCD coupling constant 
at finite temperature. We use the following 
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notations henceforth, 



,LO 



NLO 



N N f 



'D 



NP 



1 D 



,LO 



"D L ° + c N g'T + d N , Nf g i T 



Ng 2 T 

47T 



,LO 



ln> 



g 2 T 



lAm L D ° 



(4) 



where m£ is the Debye mass obtained by fit- 
ting the (colour-singlet) free energy in lattice 

QCD m. 

In the weak coupling (g << 1) regime, the 
soft scale {mo — gT at the leading-order) re- 
lated to the screening of electrostatic fields 
is well separated from the ultra-soft scale (~ 
g 2 T) related to the screening of magnetostatic 
fields. In such regime it appears meaningful to 
see the contribution of each terms in the the 
Debye mass (Eq. [TJ separately. But when the 
coupling becomes large enough (which is in- 
deed the case), the two scales are no longer 
well separated. So while looking for the next- 
to-leading corrections to the leading-order re- 
sult from the ultra-soft scale, it is not a wise 
idea to stop at the logarithmic term (as men- 



tioned in the notation 



,NLO 
l D 



), since it be- 



comes crucial the number multiplying the fac- 
tor 1/g to establish the correction to the LO 
result. In fact we found that the Debye mass 



in the NLO term 



.NLO 
l D 



) is always smaller 



than than the LO term (m^, ) because of the 
negative (logarithmic) contribution (log(l/g)) 
to the leading-order term, while the full cor- 
rection (all g 2 T terms) to the Debye mass re- 
sults positive. So, we will employ only three 
forms of the Debye masses viz. leading-order 
term/perturbative result (to^°), full (non- 
perturbative) corrections to the leading-order 
term (m^ p ), and lattice parametrized form 
{mjj) to study the dissociation phenomena of 
quarkonium in a hot QCD medium in this 
work. 

We now proceed to investigate in-medium 
modifications to heavy-quark potential and 
its application to determine the binding 
energy and dissociation temperature of the 
heavy-quark bound states. 



The in-medium heavy-quark po- 
tential: 

Let us now turn our attention to study the 
medium modifications to heavy quark poten- 
tial at T = which is considered as the Cornell 
potential, 



V(r) = — 
r 



(•5) 



where a and a are the phenomenological pa- 
rameters. The former accounts for the effec- 
tive coupling between the heavy quark pairs 
and the latter gives the string coupling. The 
medium modifications enters in the Fourier 
transform of the heavy quark potential as fol- 
lows: 



V(k) 



V(k) 
e(k) 



(6) 



where e(k) is the dielectric permittivity given 
in terms of the static limit of the longitudinal 
part of gluon self-energy [23[ 



1 + 5^(^20^1 



Ml 
k 2 



,(7) 



The result for the static limit of the dielectric 
permittivity is the perturbative one. If one as- 
sumes that huge non-perturbative effects (like 
the string tension) survive 1 above T c , the same 
could be true also for such a dielectric func- 
tion. So, there is a caveat that this (linear) 
relation of dielectric function (e) on M|, may 
also pick up modifications due to the presence 
of non-perturbative effects above the decon- 
fincmcnt point. To get rid of the complexity of 
the problem, we put all the non-perturbative 
effects (including the non-zero string tension) 
together in the effective charge (2<j/m 2 D ) of 
the medium modified potential which further 
depends of the Debye mass. The quantity 
V(k), the Fourier transform (FT) of the Cor- 
nell potential reads [24| : 



V(k) = -^(2M^- 



-ia 



V2irk 4 ' 



(8) 



Substituting Eq.© and Eq.© into Eq.© 
and evaluation of the inverse Fourier- 
Transform of the RHS of Eq.© one obtains 
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the r-dependence of the medium modified po- 
tential. The expression thus reads 



V(r,T) = ( 



2a 



a) 



exp (— mr>f) 



""D 

2a 



' D 



2a 



mo 



araj) 



(9) 



This potential has a long range Coulombic 
tail in addition to the standard Yukawa term. 
After taking the high temperature limit, the 
above potential takes the form: 



V(r) 



2a 



' D 



arriD 



(10) 



The above form (apart from a constant term) 
is a Coulombic type as in the hydrogen atom 
problem by identifying the fine structure con- 
stant e 2 with the effective charge 2a/m 2 D . 
Since the Debye mass mo always increases 
with the temperature, the effective charge 
2a/m 2 D gets waned as the temperature is in- 
creased. This makes the potential too shallow 
to bind Q, Q. This results the melting of the 
bound states. It is important to note here the 
difference between the screening of the charge 
and the screening of the range of the poten- 
tial. In the usual picture adopted to study the 
dissociation of quarkonia through the Cornell 
potential, linear term vanishes above the crit- 
ical temperature because string tension van- 
ishes. So, above the critical temperature, the 
only nonvanishing term in the potential is the 
attractive coulomb term which gets screened 
in a Yukawa form making the potential short- 
ranged. If the range of the potential becomes 
too short compared to the Bohr radius it will 
be dissolved into its constituents. However, 
in our case, this dissociation happens due to 
screening of the charge, not due to screening 
of the range of potential. Note that the con- 
stant terms in the potential (Eqj9]) are needed 
in computing the masses of the quarkonium 
states. It is equally important while com- 
paring our effective potential (EqJ5)) with the 
free energy in lattice studies (discussed be- 
low). However, the constant terms are not 
needed while comparing the values of the dis- 
sociation temperatures obtained in our model 



with the values in the lattice spectral studies. 
This is due to the different criterion is imposed 
to evaluate the dissociation temperatures (dis- 
cussed in the next section). 

We need to mention that our in-medium ef- 
fective potential V(r, T) in Eq.© agrees qual- 
itatively (and also quantitatively) with the 
singlet part of the free energy in the lat- 
tice QCD [25j]. This have been checked by 
plotting V(r, T) with rT for a fixed value of 
T/T c =3.32 [ill. 

We shall now systematically study the 
effects of perturbative and non-perturbative 
interactions on the binding energies and 
dissociation temperatures of quarkonium 
states in a hot QCD medium by employing 
the three form of the Debye mass(Eq[3]). 
In addition, we take advantage of all the 
available lattice data, obtained not only in 
quenched QCD (Nf = 0), but also including 
two and, more recently, three light flavors. 
We are then in a position to study also the 
flavor dependence of the dissociation process, 
a perspective not yet achieved by the parallel 
studies of the spectral functions, which are 
only available in quenched QCD. 

Binding energy and dissociation 
temperatures 

Binding energy of a quarkonium state at 
zero temperature is defined by the energy dif- 
ference between the mass of the quarkonium 
and the open charm/bottom threshold. At fi- 
nite temperature, the binding energy is de- 
fined as the distance between the peak po- 
sition and the continuum threshold, Ebi n = 
2m C;fc + Voo(T)_— M with M being the res- 
onance mass [26[. However, our definition is 
the conventional one viz. the 'ionization po- 
tential' in the atomic physics. 

Finally, Schrodingcr equation for the poten- 
tial (Eq fTU)) gives the energy eigen values for 
the ground state and excited states, viz. J/ip, 
tp' ', T, T' etc. for the charmonium and bot- 
tomonium states, by the Bohr's formula: 



&n — 9 



_ m n cr 

E i = -zhr (u) 



Thus E n = -E u -2± - • . are the 
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allowed energy levels of QQ bound states. 
These energies arc known as ionization po- 
tentials/binding energies for the n-th bound 
state. Thus, the binding energy becomes a 
temperature-dependent quantity through the 
Debye mass and it decreases with the temper- 
ature. 

However, there are other states in the char- 
monium and bottomonium spectroscopy viz 
Xc's and XbS and the binding energies for 
them are obtained from a variational treat- 
ment of the relativistic two-fermion bound- 
state system in quantum electrodynamics [27| . 

Figures 1-2 show the variation of binding en- 
ergy with temperature (in units of T/T c ) for 
the J/-0and Tstates, respectively. Different 
curves in the figure denote the choice of De- 
bye mass in Eq. (0J used to calculate the bind- 
ing energy from Eq. (fT2")l or Eq. lflU)) . We con- 
sider three cases for our analysis: pure glu- 
onic medium, 2-flavor and 3-flavor to see the 
flavor dependence of dissociation pattern in 
QCD medium. 

The binding energy of the J/ip(m 640 MeV) 
is considerably larger than the typical non- 
perturbative hadronic scale Aqcd- As a 
consequence, pcrturbative term (leading-order 
term) in the Debye mass takes care varia- 
tion of binding energy with temperature. The 
same argument holds good for T also. Once 
we switch-on the non-perturbative contribu- 
tions in the Debye mass through the coeffi- 
cients c N and d]si^N f i the Debye mass becomes 
so large that the binding energies of all the 
quarkonia becomes too small (even at the tem- 
perature 100 MeV) compared to their ground 
state binding energies. Indeed, the value of 
the Debye mass after inclusion of nonper- 
turbative terms is approximately three times 
larger than LO and NLO results and twice as 
large as lattice parametrized Dcbyc mass near 
2T C . The temperature dependence of the bind- 
ing energy for other quar konium states is stud- 
ied in length in Ref. fla ]. However, this is not 
the complete story, the situation may change 
once the 0(g 4 T) non-perturbative contribu- 



tions to Debye mass are evaluated and utilize 
to estimate the binding energy for quarkonia 
states . 

Thus the study of temperature dependence 
of binding energy will help us to determine 
the dissociation temperatures of the quarko- 
nium states in thermal medium. Mocsy and 
Petreczky (2(| have defined the dissociation 
temperature as the temperature above which 
the quarkonium spectral function shows no 
resonance-like structures, meaning that par- 
ticular state is dissolved. 

Physically (thermal) dissociation of a 
bound state in a thermal medium can be ex- 
plained as follows: when the binding energy 
of a resonance (viz. J/tp) state drops below 
the mean thermal energy of parton the state 
have become feebly bound and thermal fluc- 
tuations can destroy it by transferring energy 
and exciting the quark anti-quark pair into its 
continuum. So, if the binding energy of a cc 
or bb state at some temperature becomes equal 
or smaller than the mean thermal energy then 
the state is said to be dissociated. Since the 
(relativistic) thermal energy of the partons is 
3T hence the dissociation temperature Tjj of 
the n-th QQ bound state will be determined 
by the condition: 

~2 \ = ^ T D (13) 

The above condition gives the dissociation 
temperatures after inserting expression for the 
Debye mass displayed in Eq. ([4]) . However, the 
choice 3T is not rigid because even at low tem- 
peratures T <T C (say) the Bosc/Fcrmi distri- 
butions of partons will have a high energy tail 
with partons of mechanical energy e > \E n \. 

While determining the temperature depen- 
dence of the binding energy and dissociation 
temperatures the string tension is chosen to 
be o = 0.184GeV 2 . The dissociation temper- 
atures for the ground (IS), first excited states 
(2S), Xci an d Xb (IP) of cc and bb are listed 
in the Tables 1 and 2 for the Debye masses in 
the leading-order and the lattice parametrized 
form, respectively. We do not put up the 
list for the dissociation temperatures with the 
non-pcrturbative form of the Debye mass be- 




FIG. 1: Dependence of J/ip binding energy on temperature 



N-f=0 N_f=2 N_f=3 




FIG. 2: Dependence of T binding energy on temperature 



TABLE I: Dissociation temperatures for various 
quarkonia (in unit of T c 



r LO 

for m D . 



Quarkonium state 


Pure QCD 


N f = 2 


JV> = 3 


J/V 


1.1 


1.3 


1.2 


i>> 


0.8 


0.9 


0.9 


Xc 


0.9 


1.1 


1.0 


T 


1.4 


1.7 


1.6 




1.0 


1.2 


1.2 


Xb 


1.1 


1.3 


1.2 



TABLE II: Dissociation temperatures for various 
quarkonia (in unit of T c ) for mfj. 



Quarkonium state 


Pure QCD 


N f =2 


N f = 3 


J/V> 


0.8 


0.9 


0.9 


i>' 


0.5 


0.7 


0.6 




0.6 


0.7 


0.7 


T 


1.0 


1.2 


1.2 




0.7 


0.9 


0.8 


Xb 


0.7 


0.9 


0.9 



cause the values obtained are too small to ex- 
plain physically. We have taken the values of 
critical temperatures (T c ) 270 MeV, 203 MeV 
and 197 MeV for pure gluonic. 2-flavor and 
3-flavor QCD medium, respectively (28j . 



The dissociation temperatures obtained 
with the Debye mass in the leading-order 
(Table 1) is always larger than the disso- 
ciation temperatures with the Debye mass 
parametrized in lattice QCD (m£) (Table 
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2). for both charmonium and bottomonium 
states because of the larger value of the Dcbyc 
mass in lattice compared to LO values. The 
results shown in the Tables 1-2 lend suppor t 
with the recent lattice predictions [26|, l29l ]. 
The upper bound of the dissociation temper- 
atures could be obtained if average thermal 
energy is replaced by ~ T. 

Absorption by nucleons and co- 
movers 

So far we have discussed the fate of quarko- 
nia only when the presence of quark gluon 
plasma is considered. It is very well estab- 
lished that there are several aspects like ini- 
tial state scattering of the partons, shadowing 
of partons, absorption of the pre- resonances ( 
\QQg > states) by the nucleons before they 
evolve into physical quarkonia, and also dis- 
sociation of the resonances by the comoving 
hadrons. It has been argued that the absorp- 
tion by co-moving hadrons will be important 
for ijj' , due to its very small binding energy, 
while for more tightly bound resonances it 
may be weak. 

Let us briefly comment on them one-by-one. 
Shadowing of partons should play an impor- 
tant role in the reduced production of quarko- 
nia, especially at the LHC energies. It is clear 
that if shadowing is important, we shall wit- 
ness a larger effect on J/tp than on T, because 
of the smaller values of the x for gluons. At 
the same time, the effect of shadowing should 
be similar for different resonances of the char- 
monium (or bottomonium), as similar x values 
would be involved for them. 

The absorption of the pre-resonances by the 
nucleons is another source of pt dependence. 
It is important, to recall once again that as the 
absorption is operating on the pre-resonance, 
the effect should be identical for all the states 
of the quarkonium which are formed. 

This is a very important consideration as it 
is clear that if we look at the ratio of rates 
for different states of J/ip or the T family 
as a function of pt, then in the absence of 
QGP-effects they would be identical to what 
one would have expected in absence of nuclear 
absorption and shadowing, providing a clear 



pedestal for the observation of QGP. 

Thes another aspect of px dependence 
which needs to be commented upon. The 
(initial state) scattering of partons, before 
the gluons of the projectile and the target 
nucleons fuse to produce the QQ-pair, leads 
to an increase of the < p T > of the reso- 
nance which emerges from the collision fvj\ . 
The increase in the < p\ >, compared to 
that for pp collisions is directly related to 
number of collisions the nucleons are likely 
to undergo, before the gluonic fusion takes 
place. This leads to a rich possibility of 
relating the average transverse momentum 
of the quarkonium to the transverse energy 
deposited in the collision (which decides the 
number of participants and hence the number 
of collisions). Considering that collisions with 
large Et may have formation of QGP in 
the dense part of the overlapping region, the 
quarkonia, which are produced in the densest 
part (and hence contributing the largest 
increase in the transverse momentum) are 
also most likely to melt and disappear. This 
may lead to a characteristic saturation and 
even turn-over of the < p T > when plotted 
against Et when the QGP formation takes 
place. In absence of QGP, this curve would 
continue to rise with Et- 

Conclusions and Outlook 

We have reviewed different theo- 
ries/mechanisms of dissociation phenomena 
of heavy quarkonia in different stages of 
system produced in the relativistic nuclear 
collisions. Because the system formed just 
after the collision is not static and (ther- 
mally and chemically) equilibrated, it is 
rapidly expanding, hadronize, and finally 
produces many pions, photons, leptons etc. 
detected at the detectors. We discussed 
the dissociations mainly in three stages: 
initial state scattering/absorption, plasma 
interactions, hadronic/comover absorptions. 
However, we gave much emphasis on the 
second stage of the dissociation: dissociation 
in the deconfined medium/plasma. Apart 
from a comprehensive survey of different 
approaches in lattice QCD and potential 
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based studies, we present our very recent 
work on the dissociation of quarkonia in a hot 
QCD medium by investigating the in-medium 
modifications to heavy quark potential. This 
is something new because, in our formalism, 
medium modification results the (dynamical) 
screening of the color charge in contrast to 
the screening of the range of the potential in 
the usual screening picture. 
The screening of the effective charge, in turn, 
causes the energy of the quarkonium state 
depends on temperature. We have then 
systematically studied the temperature de- 
pendence of the binding energy of the ground 
(15) state, first excited (25) state, and IP 
states of charmonium and bottomonium in 
the pure and realistic QCD medium. We then 
determined the dissociation temperatures 
with the Debye mass in leading-order m^,° 
and in the lattice parametrized form. The 
results are reasonably close to the finding of 
other theoretical works based on potential 
models [2f|. On the other hand these values 
are significantly smaller than the predictions 
made by others [1, Hi| . 

In the end, we conclude that all the well 
explored and yet non-QGP effects need to 
be accounted for, before we can begin to see 
the suppression of the quarkonium due to 
the formation of QGP. It seems that this 
has been achieved at least at the SPS energies. 
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